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abstract 

The collapse of baryons into the center of a host dark matter halo is accompanied by radiation that may be 
detectable as compact ( <10 kpc) UV-continuum and Lya emission with Lya luminosities as high as ~ 10'*"- 10"*-' 

erg s"' in halos of mass M = 10" - lO'^M©. We show that the observed equivalent width (EW) of the Lya line 
emitted by these cooling clouds is EW <400 A (restframe). These luminosities and EWs are comparable to those 
detected in narrowband surveys for redshifted Lya emission. The rest-frame ultraviolet of Lya emitting cooling 
clouds radiation may be dominated by two-photon transitions from 2s^ls. The resulting spectrum can distinguish 
cooling clouds from a broad class of young star forming galaxies. 

Subject headings: cosmology: theory-galaxies: high redshift-galaxies: formation-radiation mechanisms: 
general 
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1. INTRODUCTION 

Gas collapse into the center of a host dark matter halo re- 
quires the release of gravitational binding energy in the form 
of cooling radiation. The efficiency of gas cooling therefore 
plays an important role in the formation of stars and galax- 
ies. For example, it has been shown that trapping of cool- 
ing ra diation can affect the fragmentation of a collapsing gas 
cloud JSpaans & Silkll2006l) . Understanding cooling radiation 
is therefore essential in our understanding of galaxy formation. 
A lar ge fraction of c ooling radi ation niay ern erge as Lya emis- 
sion dHaiman et al.ll20 00: Fardal et al.l 200 lb. possi blv accom- 
panied by emission from helium (lYang et al.ll2006:l . Galaxies 
may therefore be visible as extended Lya sources before star 
formation occurs. This realization coincided with the discov- 
ery of two large Lya 'blobs' (LABs) bv Steidel et al.1 fcOOOP l 
Naturally, LABs have been associate d with cooling radiat ion. 

Since the discovery of LABs by ISteidel e t alT ('2000') sev- 
eral t e ns of new LAB s have been discovered (Matsuda et al. 
120041: iDev et"an 120051: ISaito et al.l 120061: lOuchi et al.l l2008bl) . 
Various gr oups have investigated what powers Ly a emission in 
LABs (e.g. lMatsuda et al.l2006l:ISaito et al.120071) . Recent stud- 
ies suggest that the majority is likely powered by (obscured) 
stars a nd/or AGN (e.g. Basu-Zvc h & Scharf 2004; Matsuda etall 
2007 ), or by a superwind (e.g. Taniguchi & Shioval2000l : lMori et al.l 
20041) . However, in some LABs it is possible to rule out all 



these mechanisms, as each is expected to produce characteris- 
tic metal lines at flux levels that are not detected. This has led 
to the conclusion that cooling ra diation may ha ve been detec ted 
dNilsson et al.ll2006l:ISaito etall 2007: Smith & Jarvisll2007"l) . 

Cooling radiation may also manifest itself as compact Lya 
emission (rather than spatially extended) with luminosities com- 
para ble to those o bserved in high-redshift Lya emitting galax- 
ies. iHaiman et al. (2000) pointed out that ~ 30% of their to- 
tal computed coohng luminosity may originate from the cen- 
tral <10% in radius. Furthermore, iFardal et all ( 1200 lb and 
lYang et aLJ ((2006) found that spatially extended cooling gas clouds 
in 3-D simulations contain dense compact clumps (< 10 kpc, 
«H,max ^ 50 cm"-') which dominate the total Lya cooling emis- 

' Extended Lyo emission is quite common around radio galaxies, and is 
thought to be powered by some sort ofjet-IGM interaction (e.g. Chambers et al. 
.1990; .Villar-Marti'n .2007.) . However, the LABs discovered bv .Steidel et al, 
120001) have no associated radio-galaxies. 



sion. Since it is plausible that star formation takes place in- 
side these densest clumps, a major challenge lies in determin- 
ing which process provides the biggest contribution to the Lya 
luminosity emerg ing from clumps: cooling or young stars (e.g. 
iKeres et al]|2005L for a more detailed discussion). 

iDiikstra et al.l (12006*) quantified the properties of cooling ra- 
diation by performing Monte-Carlo Lya radiative transfer cal- 
culations in collapsing optically thick gas clouds, and found 
that Lya cooling radiation should be systematically blueshifted. 
In practice it may prove difficult to actually observe this blueshift, 
as this requires the detection of other lines, which are often ab- 
sent. Dijkstra et al. (2006) also showed that Lya cooling radi- 
ation may reveal itself through a frequency dependent surface 
brightness profile: the surface brightness profile of the reddest 
Lya photons is expected to be flatter than that of the bluest Lya 
photons. However, the strength of this signal is model depen- 
dent and may be too subtle to be detectable. Therefore, it is 
desirable to have more diagnostics that characterize cooling ra- 
diation. 

The main goal of this paper is to provide new diagnostics to 
identify (compact) cooling clouds in data generated by existing 
and upcoming telescopefl We compute the continuum emis- 
sion of cooling clouds redward of Lya, and show that this is 
dominated by two-photon transitions from 2s^ls. The spec- 
trum of this continuum has a peculiar shape which may distin- 
guish cooling clouds from young star forming galaxies. 

The outline of this paper is as follows. In § |2]we illustrate 
that cooling radiation may be detectable as compact Lya emis- 
sion. This is followed by a calculation of the spectrum of cool- 
ing radiation at wavelengths redward of Lya in § [3] Finally, in 
§ |4] we discuss the results of this paper, and present our main 
conclusions. 

^ The term 'Ly« coohng radiation' that is used throughout this paper refers 
to both real (i.e. collisionally excited) cooling radiation, as well as Lya ra- 
diation that is emitted as a by-product of recombination cooling (even though 
Lya radiation that is produced in radiative cascades following recombination is 
technically not cooling radiation, since it does not transport any heat from the 
gas). In both cases, the Lya radiation is emitted in a 'proto-galactic' gas cloud 
that is still collapsing to from stars. Furthermore, for neutral fractions as low 
as ,rHi = 10^, the 'real' Lya coohng rate due to collisional excitation, ALya 
is equal to the Lya emission rate due to recombination (for Tgas = 3 X 10^ K, 
see Appendix|B). Hence, even recombining H II regions can emit a significant 
amount of Lya cooling radiation. 
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Fig. 1. — The collapse of a Mto( = 10 Mq halo in a 1-D hydrodynamical simulation and its accompanying Lyo cooling radiation ai'e shown. The top panel 
shows the radial evolution of mass shells initially enclosing logMtoi/MQ = 10.5, 1 1 and 1 1 .5 (see labels). The lower horizontal axis shows the age of the Universe in 
Gyr, while the upper shows redshift. The lower three panels show the total Lyo luminosity outside a sphere of radius R as a function of R (see text for a justification 
of this pecuhar choice) at three times. The total Lya luminosity derived from the 1-D code is likely underestimated (see text). For comparison, we also show the 
luminosity in recombination (dotted line) and bremsstrahlung {dashed line) cooling radiation. These processes are clearly subdominant. In the early stages of the 
collapse {lower left panel) shells at all radii contribute to the total Lyo luminosity. However, during later stages in the collapse the Lya luminosity is dominated by 
emission from the central few kpc. This Figure illusti'ates a property of cooling radiation that is not unique to this spherical collapse model: namely that the total 
Lya luminosity of a cooling cloud may be dominated by emission from the central few kpc. 



2. COMPACT LYa COOLING RADIATION 

We use a modified version of a one-dimensional (spheri- 
cally symmetric) h ydrody namical code originally written by 
iThoul & Weinberg d 19951) . For a detailed descr i ption of the 
code the reader is referred to iThoul & Weinberg (|l995l) . The 
modifications are described in iDijkstra et alj (|2004|) . The goal 
of this section is to emphasis the possibility that the Lya lumi- 
nosity of compact cooling clouds is comparable to that of star 
forming galaxies. 

We focus on one particular simulation in which a dark matter 
halo of mass logMtot/M0 = 11 .5 collapses at Zc = 4.5. Such an 
object is expected to have a Lya luminosity of L^ ~ a few 10^^ 
erg s"' (Haiman et al. 2000 ). which is within the ra nge of lumi- 
nositie s that is probed by existing Lya surveys (e.g. lOuchi et al.l 
l2008ah . At redshift z = 4.5 the Universe was ^ 1 .4 Gyr old. We 
take snap shots of the simulation at intervals of 75 Myr and read 
out quantities such as gas density, velocity, cooling rates etc. 

The top panel of Figure [T] shows the radial evolution of mass 
shells initially enclosing logMtot/M0 = 10.5, 11 and 11.5 (see 
labels). The lower horizontal axis shows the age of the Universe 
in Gyr, while the upper shows redshift. Each shell initially ex- 
pands with the Hubble flow and is subsequently decelerated by 
the enclosed mass until it collapses. The lower three panels 
show the total Lya luminosity outside a sphere of radius /? as a 
function of /?, i.e. LiydR) = A-ir J^^ r^drja{r). The reason for 
plotting this quantity is that the Lya volume emissivity, ja{r), 
increases rapidly toward smaller radii. This causes the total 
Lya luminosity emitted by material enclosed within radius R 
(which is 47r J r^drjair)) to depend strongly on the choice 
of Ttnin. Furthermore, jaif) is highly uncertain at small radii, 
because various effects that are not included in the simulation 
(e.g. star formation) can affect jair), r^^^ and the gas' central 



density and velocity field. The Lya luminosity plotted in Fig- 
ure[T]is only affected at the smallest radii by these uncertainties. 

The lower left panel of Figure [1] shows that all radii con- 
tribute to the total Lya luminosity during the early stages of the 
collapse. However, at later times the density in the inner few 
kpc builds up (e.g. at time=l.l Gyr, the number density of hy- 
drogen nuclei is >0.01 cm"^ for r <4 kpc, see Appendix lAb 
which boosts the collisional excitation rate of the Lya transi- 
tion. This causes the central few kpc of the collapsing gas to 
dominate the total Lya cooling luminosity. 

Of course, the results presented above are suspicious as they 
were obtained from a spherically symmetric simulation. In a 
scenario of hierarchical structure formation halos are build up 
of mergers of smaller halos that collapsed at earlier times. Gas 
accretion is therefore likely to occur in dense clumps, which 
affects the results we obtained above in several ways: (i) Col- 
lisions between gas clumps convert the clumps' kinetic energy 
(due to infall) into thermal energy, which is then radiated away 
as cooling radiation; (ii) Clumping can boost the the collisional 
excitation rate of atomic hydrogen and hence the Lya emissiv- 
ity. 

Only (i) affects the predicted surface brightness profile: In 
the 1-D calculation, gravitational binding energy is mainly con- 
verted into kinetic energy of the infalling gas. This kinetic en- 
ergy must be (p artly) released as cooling rad iation in the ce nter 
of the halo (see lBirnboim & D ekel 2003; K eres et al.ll2005L for 
more detailed discussions on this topic). Therefore, most of the 
gravitational binding energy can only be released in the center, 
which causes cooling radiation emerging from a 1 -D simulation 
to peak too much toward the center. However, iHaiman et al.l 
(2000) showed that clumpy accretion also results in surface 
brightness profile that is peaked towards the center, and that 



the inner 10% of the gas (in radiu s) emits ^ 30% of to tal Lya 
luminosity. In the calculations of iHaiman et alj ( l2000l) gravi- 
tational binding energy is emitted continuously while the gas 
streams tow ards the cen t er of t he halo. Furthermore, 3-D sim- 
ulations by iKeres et all (l2005h have shown that the gas' in- 
fall speed increases down to ^ 0.5 virial radii after which th e 
infall velocity slowly decreases (also see I Wise & Abell [20071) . 
In this case gravitational potential energy is radiated away in 
a mor e compact region than that described by Haiman et al.l 
( I2OOOI) . Therefore, we conclude that the actual compactness of 
Lya cooling radiation is probably bracketed by the 1-D sim- 
ulation shown abov e and the analytic estimate presented by 
IHaiman et all (l200tf). 

The lack of clumps causes the total Lya luminosit ies shown 
in Figure[r|to be lower than the estimates presented bv lHaiman et al.l 
(I2OOO), which only depended only on the energy budget of the 
gas. Haiman et al. (2000) found that the total Lya luminosity 
due to cooling can be as large as 10"^^- 10"^"* erg s"' in halos in 
the mass range logMtot/M0 = 11-12. The compactness of the 
emission implies that a cooling blob <1 arcsec in diameter may 

emit Lya luminosities exceeding lO'^"' erg s"'. This Lya lumi- 
nosity is compara ble to that observed hig h-redshift Lya emit- 
ting galaxies (e.g. lKashikawa et alj|2006l ). It is worth pointing 
out that resonant scattering of Lya is not expected to flatten the 
observed surface brightness profile significantly (Dijkstra et al. 
2006, hence these source also appear compact on the sky). 

It should be pointed out that in the hierarchical structure for- 
mation model, the progenitor halos that merge into the halo of 
interest, are likely substantially more massive that the minimum 
halo mass that is capable of cooling gas by collisionally excit- 
ing atomic transitions (i.e. halos with 7^1, > 10"* K which cor- 
responds to Miot >1O'^M0). Hence, gas inside these halos were 
likely capable of collapsing and forming stars. It is therefore 
plausible that previously formed stars already exist in the newly 
formed halo in which gas still needs to cool and collapse to its 
center. These older stars however are not expected to contribute 
to the rest-frame UV (and Lya) luminosity of the cloud. 

3. COOLING RADIATION REDWARD OF LYa 

In the previous section we concluded that the Lya luminosity 
associated with compact cooling gas clouds can be comparable 
to that of star forming galaxies. In this section we calculate the 
spectral energy distribution of cooling radiation at wavelengths 
other than Lya, and investigate whether this can distinguish a 
cooling cloud from a star forming galaxy. 

3.1. The 2-"f Continuum: Luminosity & Spectrum 

lYang et al] (l2006h found that gas in their simulations becomes 
self-shielding for densities exceeding n >200n, where h is the 
mean cosmic number density of hydrogen (see their Fig 2). In 
the absence of a photoionizing background, coUisional excita- 
tion of hydrogen dominates the total cooli ng rate for gas tem- 
peratu res of '^ 1 -3 X 10"* K (see e.g Fig 1 of lThoul & Weinberg 
1 1995b . For higher temperatures the residual neutral hydrogen 
fraction becomes too small (xm <1 x lO""' assuming collisional 
ionization equilibrium, see e.g. House 1964, Hui & Gnedin 
1997), and this cooling mechanism is shut off. Even in the 1- 
D simulation used in § |2] which did not include self-shielding, 
collisional excitation of hydrogen dominates the cooling rate at 
radii <10 kpc (see FiglHi. 

Collisions between electrons and neutral hydrogen atoms mainly 
excite the n = 2 states (collisional excitation of the n = 3 states 
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Fig. 2. — Spectral energy distribution of cooling radiation emitted by gas 
at r ~ 2 X 10 K is shown. A prominent Lycn line (see Fig|3] the width and 
height of the line shown here are schematic.) is accompanied by a continuum 
redward of Lya. This continuum is dominated by two-photon transitions from 
2s^ls following collisional excitation of the 2s state of atomic hydrogen. For 
comparison, an uni'eddened starburst is expected to emit a spectrum of the form 
fv =constant. The units of /^ are expressed in terms of the observed Lyo line 
flux (see text). 



occurs >10 times less frequently at the temperatures consid- 
ered here). The collisional excitation rates from ls^2p and 
ls^2s are comparable as the effective collision strengths are 
comparable, il(ls,2s)/i7(ls,2p) = 0.6-0.5 for Tgas = 1-2 x 
10'* K (lOsterbrock & Ferlandll2006l Table 3.16, the tempera- 
ture dependence of this ratio is weakfl Transitions from 2s^ 1 s 
may only take place when tw o photons are emitted that have a 
combined energy of 10.2 eV (Bre it & Teller|[l940h . This two- 
photon process results in the following spectrum 

r2(is 2s) . 



17(ls,2p) 



(1) 



where M.ya is the rate at which Lya photons are emitted, 4)(v/va) 
dv/va is the probability that 1 photon is e mitted in the fre- 
quency range v±dv/2, which can be found in lSpitzer & GreensteinI 
( II95II) . The spectral energy distribution according to Eq[T]is 
shown in Figure|2] The units of /,y are chosen such that J dvf2.y{v) = 
[J7(1s,2s)/57(1s,2p)]/q,/7^. Here, /„ is the observed Lya flux 
and, and 7^ is the fraction of Lya photons that is transmitted to 
the observer through the IGM (see below). For example, Fig- 
ure|2]shows that the peak flux density in the two-photon contin- 
uum fif.m-dx = 2 x IQ~^^ f a. niO -2 /Ta) erg s"' cm"^ Hz"', where 
fa = fa,u X 10"'^ erg s"' cm"^. Also shown is a schematic 
representation of the Lya line (its shape is arbitrary). 

Other process which contribute to wavelengths redward of 
Lya include recombination and bremsstrahlung cooling. As al- 
ready mentioned however, in cooling gas at T <3 x 10"* K both 
of these processes are negligible compared to collisional exci- 
tation of hydrogen (see § 13.21 for a more detailed discussion), 
and their contribution to the continuum spectrum can be safely 
ignored. 

The prominence of the Lya line relative to its continuum is 



^ Collisional deexcitation of the 2s state suppresses the 2—7 emissivity by a fac- 
tor of (l+np/«r)"', where Mr = 1.5 X Itf* cm"^^ is the critical density of the 2s — > 
Is tra nsition , and n,, is the number density of protons ( Osterbrock & Fe rlancj 
120061) . Therefore, only gas with densities exceeding nc is expected to emit the 
two-photon continuum at levels less than that calculated above. These high 
densities are never reached in the simulation (see Appendix |A) and collisional 
deexcitation is not important in the context of this paper. 
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Fig. 3. — The emitted rest frame Lya equivalent width (EW) of cooling 
radiation as a function of the wavelength (rest frame) at which the continuum 
is measured. The EW diverges when A —> 1216 A, as the continuum flux 
density approaches zero. Note that the observed rest frame EW is likely to be 
a factor of ~ 5 — 10 lower because of scattering of Lyo photons in the IGM 
(see text). 



quantified by the equivalent width (EW) which is defined as 
EW= La/ fx- Because of the strong wavelength dependence 
of the continuum between A = 1216- 1600 A, the precise EW 
depends on the wavelength at which the continuum flux den- 
sity is measured. This is shown in Figure [3] which shows the 
emitted rest frame Lya EW of cooling radiation as a function 
of the wavelength (rest frame) at which it is measured. The EW 
diverges when A -^ 1216 A, as the continuum flux density ap- 
proaches zero. For A = 1300- 1800 A (rest-frame) the emitted 
EW is ^ 1000- 1300 A. These EWs are comparable to those 
emitted by young metal-poor galaxies A (Sc haereill2003l) . 

The observed rest-frame EW is less than that calculated above 
as residual neutral hydrogen in the intergalactic medium (IGM) 
is expected to scatter a signifi cant fractio n of the emitted Lya 
flux out of our line of sight. Dii kstra et al.l (200 7) calculated that 
only Ta ~ 0.1-0.3 of all Lya that is emitted by high-redshift 
galaxies is actually detected (7^ = 0.2 was quoted earlier since 
it corresponds to the median of this range). Moreover, Lya 
that emerges from the central collapsing clump is likely to be 
somewhat bl ueshifted (§[B, whic h results in a preferred lower 
value of Ta JDiikstra et aiul2007l) . We conclude that the IGM 
transmission is very likely 7^ <0.3, and that consequently the 
observed rest-frame Lya EW of cooling radiation is EW <400 

A. 

The observed rest-frame EW among high redshift Lya emit- 
ting galaxies can be EW ^ 100-500 A (e.g. lDawson et alj2004l: 
IShimasaku et al.ll2006l) . The existence of these large EW Lya 
emitting galaxies has led to suggestions that primordial star 
formation may have been obs erved (e.g. Malhotra & Rhoads 
l2002t iDiikstra & Wvithdl2007l) . Indeed, the results presented 
above suggest that some of these galaxies may be powered by 
cooling radiation. Figure |2] shows that the spectrum of cool- 
ing radiation redward of Lya peaks around A ^ 1600A, with 
a sharp drop down to A = 1216A. For comparison, a power- 
law spectrum of the form f\ oc A~^ (or /^ =constant) is ex- 
pected from an unreddened starburst and consistent with obser- 
vations of 27 candidate z = 6 galaxies in the Hubble Ultra Deep 
Field JStanwav et alj|2005l) . If detected, its continuum spec- 



trum could therefore distinguish cooling clouds from young star 
forming galaxies. This prospect is explored in more detail be- 
low. 

3.2. Prospects for Using 2-7 Emission as a Tracer of Cooling 

Clouds 

3.2.1. Cooling Clouds and the 2-7 Continuum 

The fate of the gas as it ap proaches the cent er of the dark 
matter halo is not known (see lKeres et al.ll2005L for more dis- 
cussion on this). The gas may be slowed down gradually, or 
may smoothly join a rotating gaseous disk that formed in the 
center of the dark matter halo. In both these scenarios the gas 
would remain cold and would emit the spectrum shown in Fig- 
ure |2] 

However, this is not the case when gas approaches r = with 
velocities of order ^100 km s"' (as was the case in the 1-D 
simulation used in §|2]i. When the gas crashes into the center, 
it is shock-heated to T ^ 10^ K and emits bremsstrahlung up 
to energies of ^ 100 eV. This ionizing radiation creates an H II 
region of < a few kpc because of the large gas densities in 
the center of the halo dBirnboim & Dekel 2003). This ionizing 
radiation could furthermore doubly ionize helium. 

This recombining gas may emit a spectrum that is very simi- 
lar to that shown in Figure|2l In Figure|4]the solid line shows the 
spectrum emitted by gas assuming case-A recombination and a 
gas temperature of T= 2 x lO'* K. In this case the bound-free 
(free-free ) emissivity is given by e = 5.3 (2.6) x 10"^^ erg s"' 
cm^ (e.g. i Qsterbrock & Ferlandll2006l) . Furthermore, the total 
emissivit y in the two-photon continuu m is e '^ 5.0 x 10"^^ erg 
s"' cm-* OOsterbrock & Ferlandll200 6l Table 4.1, and we used 
that £27 ~ 0.5eLya)- It is noted that when helium is doubly ion- 
ized inside the H II region, and/or for higher gas temperatures, 
more continuum emission is emitted at A < 2000 A (as well as 
a He II Ha line at A = 1640 A), which would weaken the appear- 
ance of the dip. The dotted line shows the spectrum resulting 
from bound-bound transitions only (which corresponds exactly 
to the case shown in Fig|2]). Technically the Lya -1- 2-7 radia- 
tion is not cooling radiation itself, but is emitted by the cooling 
cloud as a by-product of recombination cooling radiation (see 
§[Tll. Clearly, recombining gas at T <2 x 10** K emits a spec- 
trum of similar to that of 'pure' cooling radiation, especially at 
wavelengths where the spectrum differs from that emitted by 
star forming galaxies, A = 1200-2000 A. 

3.2.2. Star Forming Galaxies and the 2-7 Continuum 

Of course, recombining H II regions (of primordial compo- 
sition) surrounding stars or quasars emit the same spectrum as 
that shown in Fig H] However, in observed galaxies this nebu- 
lar emission is subdominant to the continuum radiation from 
stars/quasars, which does not exhibit the dip shown in Fig- 
ure m A notable exception is discussed by Schaerer (2002), 
who shows that nebular emission dominates the stellar contin- 
uum at wavelengths redward of Lya during the first 1-2 Myr 
in metal-free galaxies that are forming stars according to top- 
heavy initial mass function (IMF). Hence, these primordial star 
forming galaxies also exhibit the two-photon dip in their spectra 
during the first 1-2 Myr of their lifetime (see Fig 5 of Schaerer 
2002). On the other hand, these galaxies also emit He III recom- 
bination lines such as the He 1640 line, which is not the case 
for cooling radiation when the gas does not get shock heated to 
T >\Qr' K (as was discussed above in § 3.2.1). 
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Fig. 4. — The spectrum emitted by recombining gas at T= 2 X lO** K is 
shown (solid line). This spectrum includes the contribution from free-free, 
free-bound and bound-bound transitions. The dotted line shows the contribu- 
tion from bound-bound transitions only, and this spectrum is identical to that 
of cooling gas (Fig|2). Clearly, the spectra of recombining and cooling gas are 
very similar at A <2000 A. 



3.2.3. The Impact ofLya Radiative Transfer 

It is worth pointing out is that for sufficiently large columns 
of neutral hydrogen gas (A^hi >10^^ cm"^), Lya photons scat- 
ter so frequently that collisional deexcitation from the 2p level 
may become important, and the Lya cooling channel may be 
suppressed strongly. This process has been studied in more de- 
tail by Spaans & Silk (2006). However, this mechanism only 
becomes efficient when densities exceed >10^ cm"^^, which is 
much larger than the highest densities in the simulation shown 
in Figure[T](see Fig|5]in Appendix|5]). Furthermore, this effect 
enhances the expected 2-j continuum relative to the Lya line, 
and if Lya were ever detected from these ultracompact cooling 
cores, then the associated 2-7 flux is relatively brighter and thus 
easier to detect. 

So far, the role of dust in the Lya radiative transfer has not 
been discussed. One might expect little dust to exist in gas that 
is still cooling to form stars. However, as mentioned in § |2] 
in models of hierarchical structure formation, star formation is 
expected to have occurred in the less massive progenitor halos 
which merged into the halo of interest. Thus, dust is expected to 
exist in the cooling cloud. Dust would redden and suppress the 
UV-continuum, b ut to quantify this effec t theoretically is diffi- 
cult. Empirically, iBouwens et alj (l2006h found the rest-frame 
UV continuum of z = 6 Lyman Break Galaxies (LBGs) to be 
suppressed by dust by a factor of only ^ 1 .4. It is reasonably 
to expect the progenitor halos not to contain more dust than 
these more massive -and hence probably more evolved- LBGs, 
and we therefore do not expect dust attenuation of continuum 
cooling radiation to be significanQ. 

The impact of dust on Lya is less clear. In principle, reso- 
nant scattering enhances the total distance Lya photons travel 
through a medium, which makes them more prone to destruc- 
tion by dust. In this case, the 2-7 continuum would be enhanced 
relative to Lya and would therefore be easier to detect. On the 
other hand, Lya may preferentially avoid destruction by dust 

"* Gnedin et al. (2008) found the mean escape fraction of ionizing radiation in 
(adaptive mesh refinement) hydrodynamical simulations of high-redshift galax- 
ies to be (/esc) '^1%. However, the opacity to ionizing photons in these simu- 
lations was found to be dominated by neutral atomic hydrogen and helium. 



in a two-phase medium in which dust resides in cold ( T r^ 10"* 
K) neutral denser clumps e mbedded in a hot medium dNeufeldl 
II99II; [Hansen &Ohll2006 1. In this scenario, the 2-7 contin- 
uum would be attenuated (by a factor of < 1 .4 as we argued 
above), but the Lya would be freely transmitted. It is therefore 
good to keep in mind that our quoted Lya luminosities have 
not corrected for dust attenuation, but that this correction is not 
necessarily needed. 

4. DISCUSSION & CONCLUSIONS 

The collapse and condensation of baryons into the center of 
their host dark matter halos is accompanied by emission of ra- 
diation that may be detectable as compact UV-continuum and 
Lya emission with Lya luminosities as high as ~ 10"*^- 10"*-^ 
erg s~' in halos of mass M = 10^' - IO'^Mq. These luminosi- 
ties are comparable to those detected in narrowband surveys for 
redshifted Lya emission. 

Cooling clouds emit continuum radiation redward of Lya 
that is dominated by two-photon transitions from 2s^ls. The 
emitted Lya equivalent width (EW) relative to this continuum 
is EW^ 1000- 1300 A (restframe), which is comparable to the 
EW emitted by galaxies forming either metal-poor or very mas- 
sive stars. The IGM is likely to transmit only <30% of this Lya 
flux and the observed rest frame equivalent width is therefore 
EW <400 A. This is comparable to the observed equivalent 
widths among high-redshift Lya emitting galaxies. It is there- 
fore possible that existing surveys have detected compact Lya 
emission that was emitted by cooling clouds. 

Because of the dominance of the two-photon continuum, cool- 
ing clouds emit a spectrum that peaks at A '^ 1600A (rest-frame), 
and sharply drops down to A = 1216A. This peculiar spectral 
shape distinguishes cooling clouds from observed young star 
forming galaxies, whos e continuum is dominated by radiation 
from stars/quasars (e.g. lStanwav et al] |2005-). However, this is 
not true for metal-free galaxies that are forming stars according 
to top-heavy initial mass function. In these primordial galax- 
ies, nebular emission may dominate the continuum redward of 
Lya during the first few Myr of their lifetime (Schaerer 2002). 
On the other hand, these galaxies also emit He III recombina- 
tion lines such as the He 1640 line, which is not the case for 
cooling radiation when the gas is not shock heated to T > 10"* K 
(see § 13.21 ). Furthermore, deeper Lya observations of compact 
cooling clouds should find weaker emission on larger scales. 
However, the same applies to Lya emitting star forming galax- 
ies, which are likely surrounded by faint Lya halos due to res- 
onant scattering in the surrou nding IGM (e.g. iLoeb & Rvbickil 
119991: iDiikstra & Loebll2008l) . This scattered Lya emission is 
not accompanied by continuum radiation redward of the Lya 
line, which could distinguish it from spatially extended Lya 
cooling radiation (polarimetry may provide additional insights, 
see Dijkstra & Loeb 2008). 

Ongoing and future surveys for redshifted Lya emission aim 
to detect Ly a emitters out to redshifts as high as z ^ 6- 12 (e.g. 
IStark et aLi 2007: Nilsson et al. 2007, and references therein). 
These observations will provide new insights into the epoch of 
reionization and galaxy formation at high redshift. Primordial 
galaxies -like cooling clouds- may emit very large EW Lya 
emission lines (e.g. ISchaereii 12003b . In order to measure the 
EW of the Lya line and to determine whether primordial galax- 
ies are detected, the continuum redward of Lya must be de- 
tected. This is challenging; if the observed flux in the Lya 
line is fa.n (in units of 10"'^ erg s"' cm"^) and its equivalent 



width is EW (in A), then the AB-magnitude of the continuum is 
wab = 30.4 -2.5 log/a 17 +2.5 log(EW/200). For comparison, 
the Hubble Ultra Deep Field reached an 8 - cr limiting magni- 
tude in the z' filter of z^g = 28.5 (Bu nker eta l. 2004). 

Detecting both the Lya line and the continuum is not only 
desirable for identifying primordial galaxies. An additional ad- 
vantage of combining Lya and continuum measurements of 
high-redshift galaxies is that it can greatly improve the con- 
straints on the epoch o f reionization (e.g. Kashikawa et al. 2006: 
iMcOuinn et al.l2007F) over constraints obtained using Lya alone. 
The work presented in this paper suggests that if primordial 
galaxies can be detected and identified, then the same applies 
to cooling clouds. Moreover, pri mordial galaxies may only 
be luminous in Lya for <10 Myr fJMalhotra & Rhoads[|2002l: 
iDiikstra & Wvithell2007l) . Coohng clouds could be luminous 
in Lya for longer than 10 Myr, which would make them more 
easily detectable than primordial galaxies. The discovery of 
two-photon dominated Lya emitters would offer us an exciting 
glimpse of either the process of gas cooling in newly forming 
galaxies, or of star formation in the first primordial galaxies. 
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APPENDIX A 
OUTPUT FROM THE SIMULATION 

Two snapshots of the simulation used in §|2]are shown in Fig- 
ure|5] The radial dependence of the number density of H-nuclei 
(nn, upper panel), the fraction of neutral atomic hydrogen gas 
(xhi = hhi/wh, central panel), and the gas temperature (T, lower 
panels) are shown at time=0.9 Gyr {solid lines) and time=l.l 
Gyr {dashed lines). Snapshots of the simulation at these times 
were also shown in Figure [1] 

Figure |5] shows that the gas temperature in the inner 10 kpc 
is 7"gas <3 X lO'* K at both times, which corresponds to the tem- 
perature regime where collisional excitation of hydrogen domi- 
nates gas cooling. Note that the snapshot at time= 1 . 1 Gyr shows 
some strange behaviour at r <2 kpc. This is related to the for- 
mation of a shock that is produced by the gas that is radially 
moving inwards at a supersonic speed, which causes the gas to 
shock-heat to T > 10^ K. As was mentioned in §|2] the simula- 
tion results are suspicious in these innermost regions. 

Lastly, we point out that the gas density is well below the 
critical density, n^ = 1.5 x lO'* cm""*, at which collisional deex- 
citation of the 2s-state becomes important. 



APPENDIX B 

LYa RECOMBINATION VS. LYa COOLING RADIATION 

Here we compare the Lya cooling rate due to collisional ex- 
citation, Alyq, with the Lya emission rate due to recombina- 
tion, R-Lya ■ Their ratio is given by 



U- 



A 



Lya 



_ 7.3 X 10-1 Vexp(=iiMX)) 



(Bl) 



T^Lya 0.68arec,B X/Wa X (1.0-Xhi) 

where a,-ec,B = 2.6 x 10-'^(r/10'^ K)""-^ cm^ s"' is the case-B 
recombination coefficient. For T = 3 x 10"* K we find 7?. = 1 .0 x 
{xY^l/\Q~'^) (for jchi -C 1), and TZ increases with T . That is, even 
in a highly ionized plasma, a significant fraction of the emitted 
Lya photons are produced following collisional excitation of 
the much rarer neutral hydrogen atoms. 
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Fig. 5. — Simulation outputs at time=0.9 Gyr (solid lines) and time=l.l Gyr 
(dashed lines) are sliown. Tlie upper/middelAower panel sliows tlie radial de- 
pendence of the gas density/ neutral fraction of atomic hydrogen/temperature. 
These quantities were used to compute the cooling rates shown in Figure [T] 
The snapshot at time=l.l Gyr shows that the simulation behaves strangely in 
the inner 2 kpc (see text). 



